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Abstract. We prove sharp pointwise t -3 decay for scalar linear perturba- 
tions of a Schwarzschild black hole without symmetry assumptions on the 
data. We also consider electromagnetic and gravitational perturbations for 
which we obtain decay rates t~ 4 , and t~ 6 , respectively. We proceed by de- 
composition into angular momentum I and summation of the decay estimates 
on the Regge- Wheeler equation for fixed I. We encounter a dichotomy: the 
decay law in time is entirely determined by the asymptotic behavior of the 
Regge- Wheeler potential in the far field, whereas the growth of the constants 
in i is dictated by the behavior of the Regge- Wheeler potential in a small 
neighborhood around its maximum. In other words, the tails are controlled 
by small energies, whereas the number of angular derivatives needed on the 
data is determined by energies close to the top of the Regge- Wheeler poten- 
tial. This dichotomy corresponds to the well-known principle that for initial 
times the decay reflects the presence of complex resonances generated by the 
potential maximum, whereas for later times the tails are determined by the far 
field. However, we do not invoke complex resonances at all, but rely instead 
on semiclassical Sigal-Soffer type propagation estimates based on a Mourre 
bound near the top energy. 



1. Introduction 

The study of linear waves on fixed black hole backgrounds has a long history 
in mathematical relativity and very recently, major progress has been made on 
various aspects of the problem, see, e.g., [46], [17], [31], [32], [13], [13], [14], [4], 
[3D], [2H], 0B], @3 > 03 to name just a few of the more recent contributions. We 
refer the reader to the excellent lecture notes by Dafermos and Rodnianski [IT] for 
the necessary background and a more detailed list of references. Understanding 
the behavior of linear waves on fixed backgrounds is supposed to be a necessary 
prerequisite for the study of the stability of black hole spacetimes in full general 
relativity, one of the major open problems in the field. The goal of this paper is to 
prove point- wise in time decay estimates for linear waves on the background of a 
Schwarzschild black hole. To be precise, let 
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be the Schwarzschild metric on (t,r,0,<p) £ R x (2M, oo) x (0,7r) x (0, 27r). Intro- 
ducing the tortoise coordinate 

x = r + 2M\og(^-l) 

reduces the wave equation ^ g ip = to the form 

zp A zp zp 

(1.1) _^ + ^__ ^ + -A S2 ^ = 
where F = 4^. Our main result is as follows: 

Theorem 1.1. 77ie following decay estimates hold for solutions ip of (|1.1[) with 
data i/j[0] = (ipo,ipi): 

(1.2) iK^-i-^^iu, < W^lK^f+^s^o^Vo^Voil^ 

(1.3) ||(x)-VWIU~ < (t)- 3 !!^) 4 ^ 10 ^^,^ 10 ^,^!)^! 

where JJ stands for the angular derivative^. Here L p := LP(JH;L P (S 2 )) and (x) :— 
(l + |af) 1/2 - 

The rate t~ 3 is well-known to be optimal for radial data, i.e., vanishing angular 
momentum, see for example [16] . The same applies to the weight (.x)~ 4 . We remark 
that Tataru [IB] has recently obtained a striking result of this flavor but for essen- 
tially smooth data (he apparently needs a large number of angular derivatives). On 
the other hand, he derives his result in the greater generality of a Kerr background 
(for small parameter a) and also obtains a Huygens principle. We expect that our 
methods can be generalized to cover these as well as other scenarios, but we do not 
pursue this here. Another result in this direction, albeit for Schwarzschild de-Sitter, 
is due to Bony and Hafner [7] . By means of a resonance expansion they prove local 
exponential decay in that setting for compactly supported data. 

Let us mention two (related) extensions of Theorem 11.11 The first extension 
concerns the type of black hole perturbation we can cover. As stated above, Theo- 
rem [Ll] applies to scalar perturbations. However, one has similar statements (but 
with better decay, see below) for gravitational and electromagnetic perturbations 
of the Schwarzschild black hole which appear as a = —3 and a = 0, respectively, 
in the Regge- Wheeler potential, see (jl.6p below. In the case of a = —3 one needs 
the data to be perpendicular to the spherical harmonics Yq and {Xe,i}^ = _ 1 , and 
for g = one needs to require orthogonality of the data to Yq. These conditions 
eliminate a gauge freedom inherent in the problem (such as changing the mass or 
the charge ]3- We can cover these other values of a for two reasons: (i) the decay 
bounds in [17] apply to them, and (ii) the WKB analysis in Section [5] which is the 
only place where a plays a role in this paper, is insensitive to this modification. 

The second extension concerns faster rates of decay. In fact, Theorem 11.11 ac- 
tually gives an arbitrary rate of decay, i.e., t~ N for any iV, provided the data are 
perpendicular to the first few spherical harmonics (the exact number depending 
on N). This follows immediately by inspection of our proof, since [T7] establishes 



The notation a± stands for a±e where e > is arbitrary (the choice determines the constants 
involved). Also, instead of (ft 10 , ft 9 ) in iTOt one needs less, namely (p' T+1 ,p' T ) where a > 8 is 
arbitrary, see the proof in Section \E\ for details. 

2 From the point of view of the decay estimates in | 17| . these values need to be excluded as 
they are precisely the ones that give rise to a zero energy resonance. 
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accelerated rates as in Price's law [37], [38] for a fixed spherical harmonic. One 
formulation of this result reads as follows: 

Theorem 1.2. Suppose that ip[0] _L Yj where Yj are the spherical harmonics on S 2 
with eigenvalues less than £(£ + 1) with £ > 0. Then one has the following faster 
rates of decay for solutions ip of (jl.ip with data ip[0] = (ipo,ipi): 

(1.4) \\(x)- m m\\oo < (t)-^\\(x)"\f n+1 d x ^f n+1 ^o,n)\\i 

The implicit constant depends on £ and n, m are sufficiently large integers which 
grow linearly in I. 

The decay predicted by Price's law is t~ 2e ~ 3 but at the moment we only ob- 
tain t~ 2i ~ 2 , see [17]. In particular, for gravitational perturbations we take £ — 2 
and for electromagnetic ones 1 = 1 leading to the decay rates t~ 6 and t~ 4 , respec- 
tively, as stated in the abstract. Note that according to Price's law one should have 
t~ 7 and t~ 5 , respectively. 

1.1. Extension to more general data. As stated, Theorems 11.11 and 1 1 . 21 require 
the initial data to vanish at the bifurcation sphere x -co. This is clearly a dis- 
advantage of the result from the physical point of view since one would like to cover 
more general perturbations. However, there exists a classical construction by Kay 
and Wald 28 which enables one to overcome this restriction. In order to explain 
this clever geometric argument, we have to briefly digress into some more advanced 
aspects of the Schwarzschild geometry. As is well-known, the Schwarzschild coordi- 
nates (t, r, 9, ip) cover only a small portion of a bigger manifold which is referred to 
as maximally extended Schwarzschild or the Kruskal extension, see, e.g., |22) . [49] . 
This is shown by introducing a new coordinate system (T, R, 8, ip) which is related 
to the Schwarzschild coordinates by 

In Kruskal coordinates the Schwarzschild metric reads 
39 M 3 

g = 2±!f_(_dT 2 + dR 2 ) + r 2 (d9 2 + sin 2 6 dip 2 ) 
r 

for R > \T\ and r is now interpreted as a function of T and R. However, the 
singularity at r = 2M (which corresponds to R — \T\) has disappeared and nothing 
prevents us from allowing all values of T and R provided that R 2 — T 2 > — 1. 
This yields the celebrated Kruskal extension. A spacetime diagram of the Kruskal 
extension is depicted in Fig. [1] and the wedge S (which consists of the two shaded 
regions in Fig. [T]) represents the original Schwarzschild manifold. The Kruskal 
spacetime is globally hyperbolic and in order to describe the Kay- Wald argument, 
we consider two Cauchy surfaces So an d Si at T — and some small T > 0, 
respectively. The intersection So fl S corresponds to the initial surface t = in 
Theorems 11.11 and 11.21 Suppose now we prescribe initial data on So (sufficiently 
regular, with sufficient decay at spatial infinity but not necessarily vanishing at the 
bifurcation sphere T = R = 0) and consider the wave equation O g ip = with these 
data. We are interested in the future development in the original Schwarzschild 
wedge S. According to the domain of dependence property of the wave equation 
(see [22]) the development to the future of Si in S (the domain V in Fig. Q]) is 
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Figure 1. The Kruskal spacetime. The two shaded regions to- 
gether comprise the original exterior Schwarzschild manifold, de- 
noted by S. The solution of the wave equation O g ip = in the 
darker shaded region V is uniquely determined by Cauchy data on 

Si n<s. 

entirely determined by the values of ip and ipT on Si n S. Now we prescribe initial 
data (/, g) on Si such that 

(1) f{R, 0, tp) = -f(-R, 0, tp), g{R, 0, tp) = -g(-R, 9, if), 

(2) (f,g) coincide with (ip,ipT) on Si n S, 

(3) (/, g) are as regular as {tp, ipr) on Si. 

It is obvious from the spacetime diagram Fig. Q] that this can be done. Then we 
consider the solution ip of O g tp = with data (/, g) on Si. By the aforementioned 
domain of dependence property we have ip = tp in T>. The key observation now 
is the existence of the discrete isometry {T, R,6,ip) n- (T, —R,9,<p) which leaves 
the line R = invariant and guarantees that property (1) of the data (/, g) is 
propagated by the wave flow, i.e., we have ip(T, R, 9, ip) = —ip{T, —R, 0, tp) which in 
particular implies ip(T, 0, 6, tp) = for all T. As a consequence, by evaluating ip and 
ipT on So, we obtain new initial data on So which vanish at the bifurcation sphere 
and lead to the same solution in T> as the original data {ip, iPt)\t, ■ If the data are 
sufficiently regular, they have to vanish at least linearly in R at the bifurcation 
sphere which yields exponential decay with respect to the tortoise coordinate x as 
x — > — oo and our Theorems 11.11 and 11.21 apply. 

We remark in passing that the discrete isometry which lies at the heart of the 
Kay-Wald argument is a very fragile property which cannot be expected to hold in 
more general spacetimes. Recently, Dafermos and Rodnianski [12] devised a more 
robust method based on vector field multipliers which is capable of extending decay 
estimates up to the horizon. 

1.2. Strategy of proof of Theorems 11.11 and 11.21 The strategy for the proof 
of Theorem 11.11 is to decompose the solution into spherical harmonics and then 
to sum the resulting decay estimates. The wave equation (jl.ip at fixed angular 
momentum turns into a wave equation in 1 + 1 dimensions, namely in the time 
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variable t and the "radial" variable x. The angular derivatives in the estimates 
(|1.2j) and (|1.3[) then arise as weights given by powers of the angular momentum. 
This procedure is not expected to yield the optimal bounds as far as the number 
of angular derivatives is concerned. The specific numbers appearing in (|1.2I) above 
are a result of the Mourre estimate approach to the "top of the barrier scattering" 
which we develop in this paper. This Mourre estimate is non-classical in the sense 
that it needs to take into account that the top energy is trapping. We deal with 
the issue by means of the Heisenberg uncertainty principle (or the ground state of 
the semiclassical harmonic oscillator). The transition from our Mourre estimate 
to the decay in time is accomplished by means of Sigal-Soffer propagation theory 
going back to [43], but the implementation we follow is [26]. The further losses in 
terms of angular derivatives in (|1.3|) as compared to (|1.2|) are due to the Bernstein 
inequality and the L°° bound on the spherical harmonics. 

It is likely that Theorem 11.11 can be improved with regard to the number of an- 
gular derivatives required by a more detailed analysis of the spectral measure (at 
fixed angular momentum) for energies near the potential maximum. This would 
involve a reduction to Weber's equation and an explicit perturbative analysis of the 
Jost solutions instead of the more indirect Mourre-Sigal-Soffer method. However, 
since it would complicate this paper we have chosen not to follow that route. We 
emphasize that the number of derivatives JJ appearing in our theorems is deter- 
mined exclusively by the analysis near the maximum energy. There is a sizable 
literature on the topic of scattering near a potential maximum and, more generally, 
on scattering in the presence of trapping energies, see for example [2], [BJ, [8], [15] . 
[23] . [34], [39], and [44] an d the references in these papers. However, we find that 
the available results in that direction are either not sharp enough for our purposes, 
or actually can be sidestepped completely with the Mourre approach we follow here. 

For energies near zero, the Mourre estimates become degenerate. Therefore we 
need to rely on a WKB-type analysis of a semiclassical problem which we describe 
in detail in Section LOl (the semiclassical parameter being h = £~ x where £ is the 
angular momentum). The main issue for low lying energies is that the errors of 
the perturbative analysis of the spectral measure (and the Jost solutions) have to 
be controlled simultaneously for all energies near zero and all small h. This was 
accomplished in [10] and [9j. 

We remark that the technical part of this paper is entirely devoted to large 
angular momenta - in other words, to the analysis of the semiclassical equation. 
In fact, for angular momenta < £ < £q where £q is large we invoke the bounds 
from [17j and [16] . The constants appearing in the decay bounds in these papers 
grow rapidly in the angular momentum (in some super-exponential fashion) . This 
precludes us from summing them in £ and necessitates the separate WKB/Mourre 
analysis of this paper. However, since the latter only applies to large £ the methods 
developed in [T7] and [TB] are of crucial importance for Theorcm ll.il 

This paper is not self-contained, as it relies on the results of [10], [9], [17] and [15] . 
Needless to say, there is a long history concerning Price's law, see [46] as well as [11] . 
We refer to those papers as well as to the introduction of our earlier paper [17] for 
a detailed list of references and more background. We would also like to mention 
that this paper as well as [TU] , [H] , [H] and [TB] , are a result of those investigations 
into decay of wave equations on curved backgrounds with trapping metrics which 
began with the surface of revolution papers [40] , [41] . 
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1.3. Angular momentum decomposition. Restricting to spherical harmonics, 
the wave equation (|1.1|) takes the form 

(1.5) dli> - d 2 x i> + V^(x)ijj = dfy + Ht,^ = 

with the Regge- Wheeler potential 



\ r(x) J \ r 2 (x) r 3 (x) J 



(1.6) Vt, a (x) = \1 

for cr = l. However, as mentioned before we allow for other values of a as well with 
the physically relevant ones being a = —3,0. We will take I large and study the 
semiclassical operator 

(1.7) H = H{h):=-h 2 d 2 x + V{x-h) 

with the normalization V(x; h) := h 2 V£. a (x), and V(x max ; K) = 1 where x max is the 
location of the unique maximum of the potential. Thus h ~ i~ x as i — > oo. The 
maximum has the property that0 

(1.8) V'(x max ; K) = 0, V"(x max ; h) ~ -1 

uniformly in h and x max ~ M. For the cosine evolution one has for x' < x 



cos 



i 

cos(<A) Im 



f + (x,\;£)f-(x',X;£) 



Xd\ 



2 r°° 

(1.9) = — / cos(fi _1 ^)Im 

Tfi 2 Jo 



f + (x,E;h)f-(x',E;h) 



EdE 



.W(f + (;E;K),f-(;E;h))- 

with f± being the outgoing Jost solutions for the original operator and the 
semiclassical one, respectively. For the latter case this means that 

{-h 2 d 2 x + V(x; h))f ± (x, E; K) = E 2 f ± (x, E; K) 

f±(x,E;h) - e ± i Ex x->±oo. 

Furthermore, we write W(f, g) := fg' — f'g for the Wronskian of two functions /, g. 
The sine evolution is given by 



2 f°° 

— / sin( fr 1 tE)e(E,x,x';K)dE 
nri Jo 



(1 1Q) sin (ty/7U^)(x,x') 

with the semiclassical spectral measure 

v W{f+{-,E]h),f-{-,E;h))- 
In order to control the semiclassical evolution we distinguish energies < E < e, 
e < E < 100 and 100 < E. Here e > is some fixed small constant which does not 
depend on h. 

The regime of large energies is relatively easy, whereas the low lying energies as 
well as those near the maximum V = 1 represent the most difficult contributions to 
analyze. For small energies we follow the analysis of [TU] and [5] which was specif- 
ically developed with this application in mind. In the former paper, the challenge 
was to carry out the WKB analysis for a smooth, positive, inverse square poten- 
tial uniformly for small h and small energies < E < e. This was accomplished 



■^Throughout, a ~ 1 means that C 1 < a < C for some constant C. 
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by means of Langer's uniformizing transformation which reduces the perturbative 
analysis to an Airy equation. We note that a novel feature was the modification 
of the potential, see Section [5] In this paper we have to go beyond QJ5] since the 
Regge- Wheeler potential exhibits inverse square decay only in the far field, whereas 
it decays exponentially towards the event horizon. This is where [9| applies, which 
develops a normal form reduction for the exponentially decaying region to the left 
of the maximum of the potential. 

As mentioned above, we do not employ a uniformizing transformation for energies 
close to the maximum V = 1; this can indeed be done, and requires a perturba- 
tion theory around Weber's equatiorQ. Instead, we prove a Mourre estimate near 
the maximum. This is somewhat unusual as the maximum energy is trapping and 
therefore needs to be excluded in the classical Mourre theory, see [20], [24]. How- 
ever, a simple application of the uncertainty principle (or the ground state of the 
semiclassical harmonic oscillator) allows one to deal with this trapping case as well. 
Of crucial importance here is that the maximum of the potential is nondegenerat^l 
(in fact, V"(x max ) < 0). Once the Mourre bound is established, we employ a semi- 
classical version of the propagation estimates of Hunziker, Sigal, Soffer [26], which 
in turn go back to the work of Sigal, Soffer [43], see Section [6] below. 



1.4. Notations. In this paper we frequently employ the notation a < b (for a, b £ 
K) meaning that there exists an (absolute) constant c > such that a < cb. We 
also use a > b and write a ~ b if a < b and b < a. Furthermore, 0(f(x)) denotes a 
generic complex-valued function that satisfies \0(f(x))\ < \f(x)\ in a domain of x 
that is either stated explicitly or follows from the context. We write Or(J(x)) to 
indicate that the respective function is real-valued. The symbol ~ is reserved for 
asymptotic equality, i.e., f[x) ~ g(x) as x — > a, where /, g are two complex-valued 
functions, means that 

hm ^ = 1. 



2. LOW LYING ENERGIES AND WKB 

In this section we bound (jl.9|) for energies < E < e. Our approach is based on 
the perturbative WKB analysis of the Jost solutions which was developed in [TO] 
and [5]. More precisely, [TU] applies to the case of x > for which the potential 
decays like an inverse square. For x < the potential exhibits exponential decay 
as x — > — oo and [9] develops the methods needed for that case. We present the 
main steps of the analysis developed in these papers but omit the most involved 
technical details so as not to disrupt the flow of the argument. 



2.1. The far field. We begin with the former case, i.e., x > 0. In fact, we shall 
apply the analysis of this section to x > xq where xq < is a fixed constant. In 



We will pursue this matter elsewhere. This approach seems needed in order to prove a 
Huyghcn's principle along the lines of this paper. 

^Technically speaking, the methods of this paper apply to potentials with a unique maximum 
satisfying V^\x) = for 1 < j < k, V^ k+1 \x) ^ 0. However, the number of derivatives J7 in 
our decay estimates would then increase with k. 
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fact, any xo < is admissible, see below. Define a modified potential 

V Q (x;h) :=V(x;H) + ^(x)- 2 

and denote by x\ (E; ft) > the unique positive turning point for anjQ < E < 
i.e., the solution of Vq(x; ft) = E 2 , x > 0. 



2.1.1. Liouville- Green transform and reduction to a perturbed Airy equation. The 
analysis of [TU] was based on the following "Langer transform" , which in turn is a 
special case of the Liouville- Green transform, see [35], [33]. See Lemma 3 in [TU] 
for essentially the same statement. 

Lemma 2.1. WitM Qo ■= V - E 2 

3 



(2.1) C - C(x, E; ft) := sign(x - Xi(E; ft)) - / y/\Q (u,E;h)\du 



xi(E;h) 

C 



defines a smooth change of variables x t— > £ for all x > xq. Let q :— —Op-. Then 
1 > °> % = C = VQ, and 

-h 2 f" + (V - E 2 )f = 

transforms into 

(2.2) -h 2 w(0 = (C + h 2 V(C,E;h))w(0 

under w = %/C 7 / — Q J f ■ Here ' = ^ and 

f/ 1 -i/ \-2 -h d2 ^ 1 
V:=~ A q (X) -q ^ 

Proof. It is clear that (|2.ip defines a smooth map away from the point x = x\ (E; ft). 
Taylor expanding Qo(x, E;h) in a neighborhood of that point and using that 
V (xi(E;H)) < implies that £(x,E;K) is smooth around x — x\ as well with 
('(xi,E;h) > 0. Next, one checks that 



* = q-*f'+-hrf, w = q-*f" + ^rf 
and thus, using -ft 2 /" = {E 2 - V)f, 



dq X , .. -I „ d 2 qi 

W f ' w = q f 



-h 2 w = q-\E 2 - V)w - h 2 q-i 



- w 



q-\-Q + h 2 {x)- 2 /4:)w-h 2 q-i 



d 2 qi 
~d^ 

d 2 qi 



-w 



dC 

= Cw(0 + n 2 ( q - 1 (x)- 2 /4- q -i^-) w 

as claimed. □ 



'It was shown in |10| that the WKB approximation can only bo applied after this modification. 

This condition will be tacitly in force throughout this section. 

'We warn the reader that what was called E in 1 101 is now called E 2 . 
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The Airy equation (|2.2j) provides a convenient way of solving the matching prob- 
lenfl at the turning point £ = 0. We remark that it was assumed in [10] that the 
potential satisfies V(x) = pLx~ 2 + 0(x~ 3 ) as x — > oo. However, the methods of that 
paper equally well apply under the weaker assumption V(x) — [ix~ 2 + 0(x~ a ) as 
x — > oo where c > 2; furthermore, one needs d k O(x^' J ) — 0{x~ tr ") for all k > 0. 
This is relevant here since the Regge- Wheeler potential exhibits a ^^^-correction 
to the leading x~ 2 -decay. 

2.1.2. A basis for the perturbed Airy equation. For the Airy functions Ai, Bi appear- 
ing below we refer the reader to Chapter 11 of [35 . To the left of the turning point 
a fundamental system of (12.21) is described by the following result, see Proposition 8 
in pjj]. 

Proposition 2.2. Let Hq > be small. A fundamental system of solutions to (|2.2[) 
in the range £ < is given by 

0i (C, E, h) = Ai(r) [1 + hai (C, E, ft)] 

MC, E, h) = Bi(r)[l + ha 2 (C: E, ft)] 

with t '.= —h~^C- Here ai,a 2 are smooth, real-valued, and they satisfy the bounds, 
for all k > and j — 1, 2, and with £o := C{ X Q> E), 

sup \d k Eaj (t,E,h)\ <E- k 
to ^ Co<C<o 

2 X[-i<c<0] + ICI 2 X[Co<C<-i] 
uniformly in the parameters < h < Hq, < E < e. 

Proof. This is essentially Proposition 8 in [TU] . The two differences are (i) we work 
with E 2 instead of E (ii) the potential has a i^^-correction to the leading inverse 
square decay rather than the x _3 -correction assumed in |10) . 

As far as (i) is concerned, we note the following. Let E be as in [TU] and assume 
E 2 := E. If \d E a{E)\ < E' 1 , then b(E) := a(E) satisfies \d E b{E)\ < E" 1 by the 
chain rule. So it makes no difference whether we work with E or E. As for (ii), we 
note that the only change to the estimates in Section 3 of [10] is in Lemma 7, where 
one needs to replace (x)~ 3 f3i(x, E) with (x)~ 3 log x Pi(x, E), cf. (3.20) and (3.24) in 
that paper. However, inspection of the proof of Proposition 8 there reveals that this 
does not affect the resulting bound in any way, see the paragraph between (4.12) 
and (4.13) there. □ 

We remark that Proposition 12.21 would fail if we had defined ( in (|2.1[) with 
Q = V — E 2 instead of Qo- In the region ( > we have a basis of oscillatory 
solutions as described by the following result, see Proposition 9 in [10] . 

Proposition 2.3. Let Hq > be small. In the range £ > a basis of solutions 
to (|2.2p is given by 

MC, E\ h) = (Ai(r) + »Bi(T))[l + ft&i(C, E; ft)] 

MC, E; K) = (Ai(r) - iBi(r)) [1 + ft6 2 (C, E; ft)] 



\d k E d c a,(C,E,h)\ <E- k 



®The usual WKB machinery requires solving two matching problems, namely between the 
Airy region and the oscillatory region on the one hand, and the Airy region and the exponential 
growth/decay region on the other hand; see for example 1331 . 
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with t := — hri^ and where b\, b 2 are smooth, complex-valued, and satisfy the 
bounds for all k > 0, and j = 1, 2 

|aI^(C^^)l<C fcJ E- fc (0- f 

|3 f fe(C,£; fi)l < c fc £r fc Jr*<cr a 

uniformly in the parameters < h < Hq, < E < e, Q > 0. 

2.1.3. XTie outgoing Jost solution for the far field. We can now draw the following 
conclusions from Propositions 12.21 and 12.31 about the outgoing Jost solution /+. 
First, recall that the Airy functions satisfy the following asymptotic expansions 
with £ = 

Ai(.x) = -^tr(l + OiC 1 )), Bi(x) = -^(1 + Oir 1 )) 
27T2a;4 7r 3 a; 4 

(2.5) Ai(-x) = - Jlt-[cob(C- ~tt)(1 + 0(r 2 )) + sin(£- J^Ot 1 )], 
Bi(-z) = -J^ [ - sin(£ - i7r)(l + 0(r 2 )) + cos(e - ~k)0(C ')] 

7T 2 X 4 4 4 

as x — » oo, see 35 . Moreover, 

1 



In what follows, 



W(Ai, Bi) = AiBi'- Ai'Bi: 

7T 



ll(B;fi) 



S+{E-h):= ^/V (y;h)-E2dy 

J Xq 

T+(E;h):=Exi(E;h)- {VV {y; h) — E 2 — E) dy 

Jx 1 (E-h) 

One checks that S+{E;h) ~ |log£| as E ->• 0+, whereas T + (E;h) -> T + (0;h), 
some finite number. Moreover, 

|d|S(£; S)| + \d%T + (E; H)\ < C k E~ k 

for all k > 1. One has£3 

(2.6) Z+fo-E) = V^*fi-M^ + *>gr*(0$j(C,-E). 

This is obtained by matching the asymptotic behavior of / + with that of "02(C) as 
x — y oc and by using the relation w = f from Lemma l2.ll We refer the reader to 
[TU] for all the details. We now connect ip2 to the basis 4>j(C> E) of Proposition! 

MC,E) = c 1 (E)MC,E)+c 2 (E)MC,E) 

where 

, . W(M;E)M,E)) (FV= W(M;E),M;E)) 
cu > w(M-:E),M-,E))' 21 ] W(M;E),M;E)) 

By Proposition 12. 2[ 

W{M; E),M; E)) = -ft-t W{M, Bi) + O(hi) = -tT^-S (1 + 0(h)) 



1( ^We suppress f« as argument in most functions, even though everything here does depend 



on h. 
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where we evaluated the Wronskian on the left-hand side at C = 0. Next, by Propo- 
sitions O and El 



W(M-,E),M-,E)) = -/T*[(Ai(0) - iBi(Q))Bi'(0) 

- (Ai'(0) - iBi'(0))Bi(0) + 0(h)] 

= -ft-l[W(Ai,Bi) + 0(h)] 

(2.7) 2 
W(1> 2 (;E)M;E)) = -R"»[(Ai(0) - iBi(0))Ai'(0) 

- (Ai'(0) - iBi'(0))Ai(0) + 0(h)] 
= -h-i[iW(M, Bi) + 0(h)] 

so that 

(2.8) c 1 (E) = l + 0(h), 02(E) = -1 + 0(h) 

where the O(-) terms satisfy \d%0(h)\ < C^hE^ k . With Co = C(xo,E) we infer 
that 

Z+^o,^) = ^e'(— + i'E^-5 9 -i (Co)MCo,E) 

= v ^e l( ^ + 4)^4n-6 (? -4(C o )[ Cl (i;) ( /, 1 (C o ,^)+c2(£;)02(Co, J B)] 



/;(^o,^) = v^e i (^+i)^H-k , (^o)9-^(Co)[V'2(Co,S)- i ^(Co)V'2(Co,-B)] 
= v ^e l (^+T)i;2n-6g4( Co )[ Cl ( £ ;)0 1 (c o ,i;) + C2 (^) ( /, 2 (c o ,£;) 
- i £ (Co) (ci (£)</>! (Co , £) + c 2 (Co , )] 
where we have used that C' = (Pi see Lemma |2. II Furthermore, 

ho) = -Co 1 + icoi* Z^ *. = o(i6,i*) 

3 (Vb(xo)-£^ 2 )5 

with 0(9|;|Co|') = 0(E~ k ) for fc > 1. From Proposition 12.21 with K denoting a 
real-valued term, 

<MCo, £) = Ai(-?riCo)(l + R (fr)) 

^2(Co,£) = Bi(-n-tco)(i + P»(n)) 

M(o,e) = -/r§Ai'(-?rtCoXi + Or(»)) + QR(ft)|c |*Ai(-jHco) 

<MCo, S) - -JT'Bi'(-fi-f Co)(l + Or(S)) + B (fc)|Co|*Bi(-7rlCo) 
which implies via (|2.5p that 

MCo,E) = (4Tr)-Hh-^Co\r i e-i h ~ 1 ^(l + Om(h)) 
M(o,E) = 7r-i(^-*]Co|)-*ef R ""l<°l f (1 + R (fi)) 
^(Co,^) = ?i-i(47r)-i(a-i|Co|)^-^" 1|fo|i (l + R ^)) 
^(Co,S) = -S-t 7 r-4(n-||Co|) i el^ 1 ^l i (l + R (S)) 
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-i i i 3 

In view of these properties and using that e~ h l c °l 2 = O r (H) where d%0 M (h) 
0(E~ k h), one obtains (with C2 as above) 

U(xo,E)=c 2l q-±(Co)M(;o,E)[(l + Om(h))+i(l + 0(h))e- 2h ~ ls +] 

^ , . 1 

.f + (x , E) = c 2l q-(Co)M(o,E) [(1 + R (h)) - i(- + 0(h))e- 2h S +] 

T + (E) x , 1 1 

where 7 = 7(1?, ft) := — y^e^ s '"t) E^hTe and with 
3 



|Co| ^ = = S+= ^/V (x;h)-E2 dx 



being the action integral defined earlier. Furthermore, it follows from Proposi- 
tions [2]2] and [23] that each differentiation in E loses one power of E (in particular, 
the 0(h) terms have this property). For future reference, we remark that 

(2-10) = -^(1 + R (hj) [1 + 0(e-^ ls +)] 

where di > is a constant (depending on xq). In particular, 



Im 



= h- 1 0(e- 2h " s ^ E ^) 



2.2. Approaching the event horizon. We now deal with the potential for x < 
Xq. Here Xq < is chosen such that the Regge- Wheeler potential (setting 2M = 1 
for simplicity) can be written as 

00 

V(x;H) = ^c n -i(K)e nx 

n=l 

as a convergent series x < xq uniformly in h £ (0,Hq\. In fact, since the Lambert 
function W(z) defined via W(z)e w ^ z>> = z is analytic on \z\ < e -1 , it follows that 
one can take any xq < 0. The coefficients have expansions in powers of h and 
we normalize such that cq{K) = 1 + 0(h 2 ). One can also check that ci(0) 7^ 0. 
The goal is to control the Jost solutions f-(x,E;h) as x — > —00 uniformly for 
(E,h) £ (0,e) x (0,ft ). 

2.2.1. Transforming the problem to a compact interval. For notational convenience, 
we switch from x to — x and consider x > \xq\. The problem is then to control /+ 
for the problem 

(2.11) -h 2 f1(x,E;h) + V(x;h)f + (x,E;h)=E 2 f + (x,E;h) 

with 

00 

V(x;h)=e- x (l + J2 c n(K)e- nx ) 
n=l 

We now transform this case into a semi-classical scattering problem on a bounded 
interval (0, yo) by introducing the new independent variable y = 2e~?. Setting 
f(x) — g(y) reduces finding the outgoing Jost solution to the equation for g(y) = 

g(y,E;h), 

-h 2 [g"( y ) + y^g'iy)] + (si(y; h) %-)g{y) = 



2n 
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with the normalization g(y) ~ (y/2)~ 2l 9: as j/ — > 0+, and with 

oo 

(2-12) n(»;fe) = l + X;. 

n=l 

analytic in \y\ < y :— 2e~ x °l 2 . Finally, setting g(y;h) := y^g(y;h) yields the 
equation 

(2.13) - H 2 g"(y, E; h) + K) - + 4£ 2 ) 2 T 2 )g( 2/ , £J; R) = 
with the normalization 

(2.14) g(y,E;h)~2 2 ^ y ^ 2 ^ 

as y — > 0+. We remark that in the case 7i) = 1 the equation (|2.13l) is a modified 
Bessel equation with a basis given by the modified Bessel functions Ii l/ (h~ 1 y) and 
Ki l ,(h~ 1 y) with v = 24L It is shown in [9] by means of a suitable Liouville- Green 
transform that this basis leads to an actual basis of (|2.13l) . We begin with the 
following normal form result from |5] which is based on a Liouville- Green transform 
(the variable z below is a rescaling of y: y — az, with a := y / fi, 2 /4 + 4E 2 ). Recall 
that / — 1 means that C _1 < f < C for some constant C. 

Lemma 2.4 ( 9 ). Let Q be as above and ao > be sufficiently small. For all 

< a < a>o there exists a C°° diffeomorphism w — w(z,a) : Io(a) := (0,a~ 1 yo) — ► 
Jo(a) := (0, a~ 1 wo(a)) where yo is as above with the following properties, uniformly 
in < a < ao : 

• w (a) ~ 1 

• w'(z, a) ~ 1 for all z e Io(a) 

• \d%diw(z,a)\ < C k ,t (z) 1+k a e for all k,£> and z e I Q (a) 

Let Hi := oT x h. Then there exists a function V2(-,a;h) such that ip solves the 
rescaled form of (|2.13p . viz. 

(2.15) - h 2 2jj"(z) + (Cl(az; h) - z- 2 )ijj{z) = 

on Io(a) iff ip(w) :— (w'(z,a))iip(z) (where w = w(z, a)) solves 

(2.16) - h\Lp"{w) + (1 - w- 2 )ip(w) = hlV 2 (w)cp(w) 
on Jo(a). Furthermore, the potential Vzi- , a; K) satisfies 



\d k a diV 2 (w,a-K)\ < C k/ (w[ 



a 



for all w € Jo(&) an d k, £ > 0. 
Proof. This is done by setting 



dw z 2 — Q,{az; H) 



dz V w — 1 

More precisely, with z t being the turning point defined by zf 2 — Q(az t ; h) = 0, this 
means that 



/W pZ 
\/ \ — v~ 2 dv = I y/ f2(au; K) — u~ 2 du 



provided z > zt and 



yj — 1 + v 2 dv = / \f— tt(au; K) + u 2 du 
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provided < z < z t . Note that w — > as z — » 0. The properties of w stated above 
are now shown by calculus. The potential Va is given by 

^(W) = (W'(Z)) 2- 



4 (w'(z)) 4 2(w'(z)) 3 ' 
We refer the reader to [9] for further details. □ 

We remark that the proof also shows that w(z) = z + 0(z 2 ) and w'(z) — l + 0(z) 
as z — > 0. One now concludes the following concerning a basis of (I2.13[) . Let 
a := ^/h 2 /4 + AE 2 . Since (|2.15[) is a rescaled form of (12 . 13[) . one can now obtain a 
system of fundamental solutions to the latter equation from a perturbative analysis 
of (|2.16j) . The modified Bessel functions h v {z) and Ki V (z), which are both analytic 
on C \ (—oo, 0], give rise to a fundamental system of the homogeneous equation on 
the left-hand side of (|2.16p . In our case v = 2^. Recall the asymptotics 

(2-17) h»[z) = ^/ 2 l\A l + 0{z 2 )) z^O. 

1 (iv + 1) 

Note that for our purposes it suffices to consider real z. Moreover, Ii V {x) grows 
exponentially asi->oo, whereas Ki U (x) decays exponentially as x — > oo. 



Corollary 2.5. Let a :— ^/h 2 /4 + AE 2 with h and E > small. There exists a 
fundamental system of (j2.13[) . denoted by (go,gi), of the form 



g (y,E;h) = (w(z)/w'(z)) 1 2l 2i§ (^)(l + hc 1 (y,E;h)) 



My,E;H) = (w(z)/w'(z))h_ 2i ^(^-)(l+hc 2 (y,E;h)) 

hi 

where w(z) = w(z,E;h) is as in Lemma \2.4\ and with z — —, K\ — — . The Cj 
satisfy for all k, £ > 0, 

\d k E dl C] {y,E-K)\<C k .,a- k 

and all < y < yo. 

Proof. This follows from two facts: (i) a basis of the homogeneous equation (|2. 16|) 
is given by 

yj 

o (w,E;h) := VwI 2i E{ — ), ^(w^E;^ := VwI_ 2 E(i—) 
K h\ h hi 

and (ii): the equations for c\ t 2 are contractive; in fact, they are given by the usual 
Volterra equation involving the homogeneous basis and the potential V 2 . For ci 
one has (suppressing E and H as arguments) 

c 1 {w) = -h\ 1 j I cf)Q 2 (v)dvV 2 (u)(f>l(u)(l + hid{u))du 

JO Ju 

which implies the desired bounds on ci via Lemma 12.41 and the well-known asymp- 
totic behavior of the modified Bessel functions. For this see [9]. □ 
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2.2.2. The outgoing Jost solution towards the event horizon. From (|2.17|) . Lemma [2.4l 
and Corollary 12.51 we conclude that 



y\h(iy/2ti)- 



2'M 



y vy ' ' ' \aJ r(l-2if) V ' 

as y —> 0. In view of (|2.14[) this implies that the outgoing Jost solution is represented 
as 

rfi iv\ 

f-(x,E',K) = y/a ^_ ih y v 9i(y,E;h) 

= ^_' {aw{z)/w l {z))^I. iv {-^-){l + 0(h)) 

for all x < xq with y — az — 2e% and v = 2j-. In particular, using the standard 
asymptotic behavior of i_j„, see [5], one obtains for x — xo that 

/_(x , E; h) = 7 _v^eA( s -( E ^)+ ir -( E > fi ))(l + 0( e ^ s -))(l + R (ft)) 

/Kio^i *) = j'_^tl-^ s -^ +iT -^(l + 0(e-i s -))(l + O m (h)) 

with constants |7_| ~ \"/'_\ ~ 1 depending on iS, ft, as well as a suitable action 
S-(E; h) which is analytic for \E\ < H with S-(E; K) > for small real- valued E 1 , 
and T_(E; H) some real- valued function of real E 1 analytic on < S. We remark 
that ^7- > 0, which is most important in Section [2^41 Furthermore, each derivative 

in E costs at most a power of hr 1 . It is important that one does not lose E^ 1 as 
in the x > case, but only 0(h~ x ) as such a loss is negligible compared to the size 
of e* s -( £ *>. 

2.3. The Wronskian of the outgoing Jost solutions. From Sections 12 . 1 . 31 and 
12.2.21 it is now a simple matter to determine the Wronskian between the outgoing 
Jost solutions. 

Lemma 2.6. Define 

S(E;h) :=S+(E;h)+S-(E;h), T(E;h) :=T + (E;h) +T_(E;h) 

One has 

(2.19) W(f+(; E; &),/_(•, E- h)) = 7o Eh' 1 e h ~^ s[ - E ' h)+lT{E > h)) (I + 0(h)) 

where 70 7^ is an absolute constant, and \dgO(h)\ < Ck h~E~ k for all < E <C 1, 
< h < 1 and k>l. 

Proof. This follows from flU) and (|2T8| . □ 

Due to the growth of the action 5 , (£'; S) one can now conclude the following 
important size estimate on the Wronskian: 

\W(f+(-, E; h), /_(-, £; ft))| ~ h^Ee*^ 8 ^ > N^E) 1 ^, N := ft" 1 

for all < E < e. More precisely, one uses that S+(E; H) = — \ogE + oq + o(l) as 
E — » 0+ uniformly in small h, whereas S-(E; h) > for small £\ Then /1 := e~ Q °. 
In other words, the Wronskian blows up as E — > 0+ as a power law with large 
power since H is very small. 
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f+(x , E; h)f_{x ,E; ft) - f + {x ,E- h)f-{x ,E; ft) 
f!_(xo,E;h) /;(x ,£;ft)n-i 



2.4. The spectral measure e(E, x, x'; ft) near the maximum of the potential. 

We now derive the contribution of energies < E < e to the desired pointwise decay 
of (11.91) in time. We shall fix x — x' = xq since this case can be treated most easily 
from the previous sections; moreover, the region near the maximum of the potential 
is in some sense the most important one. The case of general x, x' is considered in 
Section 1231 First, one has 

fJP t\ y \ f-(x ,E;h)f + (x ,E;h) 

e(E; xq, Xq; n) = 1m 

= Im 

= ftlm 

where 

a(E;H) := ft" 
From H23HD one ha£3 



f-(x ,E;h) f + {x 0l E;H) 
a(E; ft) 



- a(E; h)/3{E; ft) 
1 f-(x ,E;h) 



r_{x Q ,E-ny 



P(E;h) := ft 



f'+(xo,E;h) 
f+(x a ,E; ft) 



a{E;h) :=d [l + O{e- 2h 1 S -)}(1 + R {h)) 

where do > is a constant that depends on xo,E,H with do — 1< The 0(h)- 
terms in the numerator and denominator are not necessarily the same. Similarly, 
from (|2.10p . with a constant d\ > 0, 

P(E; h) := -di[l + 0(e- 2ri " s +)](l + O m (ft)) = -d t (ft)[l + 0( e - 2ri " s +)] 

Due to the exponential decay of V as x — > — oo, the functions f-(x,E\h) and 
f'_(x,E;h) are analytic in £ in a disk |_E| < ft. In particular, a(E;h) is an- 
alytic around i? = in the same neighborhood. Moreover, due to f-(x,E;H) = 
f-(x, —E;h), one checks that Re f-(xo, E; ft) and Im /-(xq, E; ft) are even and odd 
in E, respectively. Thus, it follows that lma(E; ft) is odd in E, whereas Rea(E; ft) 
is even. Moreover, for any k, n > 0, 

(2.20) \d k E lma(E;K)\<C Kn h n 

which follows from the fact that S-(E; ft) > uniformly in — e < E < e as well as 
the differentiability properties of S-(E; ft) in E, see 0. In view of these properties, 

e(E;x ,x ;h)-"~t '"^ 

(2.21) 

Since S + (E; ft) 



ftlm 



ftlm 



1 + d x a{E\ ft) - a(E; ft)(di + ,3(i5; ft)) 



a(E; ft) 



log E as E 



l + di<x(E;h) 
►0+, 



Ob 



-2ft" 1 S+(E;ft)^ 



; -2fi- 1 S+(£;;?!,)\ _ 



= o(£; Ar ), iv = ft- 



Moreover, the imaginary part in (|2.2ip . i.e., 



Im 



a(E;h) 



l + d ia (E;h) 



llr This a is not related to a appearing in Section l2.2l 
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is an odd function in E (and analytic near E = 0) and it satisfies the bounds 

\d k E r)(E;h)\ <C k>n h n 

cf. (|2.20p . and d^r](0;h) = for even k. It is now easy to bound (11.9[) : for any 
n > 0, and all t > 0, and any < k <C ft , 



sin (H~HE) Im 



f + (xo,E;H)f-(xo,E;h) 



W{f + {x^E-K)J.{x 0l E-K)) 



Xe(E)dE 



< 



-k 



by integrating by parts (here Xe is a smooth localizer to energies E < e). In other 
words, by taking ft sufficiently small one can achieve any rate of decay. Moreover, 
we note the important property that small energies do not present any kind of 
obstruction to the problem of summing in the angular momentum £; in fact, the 
contributions of low lying energies to the decay estimates decay rapidly in I. 

2.5. The weighted L 2 bound on the spectral measure. Here we generalize 
the analysis of Section [2T4l to allow for general x,x'. More precisely, we claim the 
following result which is a routine application of the basis representations which we 
have obtained above, cf. Section 8 in |17j . 

Lemma 2.7. Let < M -C hr x . The spectral measure as defined in ([l.lip satisfies 
the bounds. 



(2.22) sup \\{x)- k -2-d%e(E,x,x';fi){x'y 

a<E<e 



for any n>0 and any < k < M. Moreover, for any choice of x,x' G R, one has 
the property 



(2.23) 

for any < j < 



lim d 2 Je(E,x,x';h) = 



£->0+ 



M 



2.6. The decay estimate for small energies. It is now a simple matter to 
establish the desired decay estimate for (|1.9p for small energies. 



Proposition 2.8. Let < M < hr 1 be given. Then for any n > 0, < k < M , 

and t > one has the bounds 



(xy 



sin (h- 1 tE)e(E, x, x'; h)xe(E) dE(x') 



2->2 



poo 

'~ / cos (h- 1 tE)e(E, x, x' ; h)xe(E) EdE(a 
Jo 



/\-fc-i- 



<C k .nH n (t)- k 
<C k ,nh n {t)- k 



uniformly in < H < ho. 

Proof. This follows directly from Lemma l2~7l by repeated integrations by parts. □ 

We remark that these estimates immediately transfer to L 1 —} L°° bounds by 
means of Bernstein's inequality. In fact, one can also establish the following re- 
sult with the optimal (x)~ k weights by means of a more careful treatment of the 
oscillatory integrals as in [17] . 
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Proposition 2.9. Let < M < hr 1 be given. Then for any n > 0, < k < M , 

and t > one has the bounds 



sin (frHE)e(E, x, x'; h) X e{E) dE(x'y k < C k , n Ti n {ty k 

1— >oo 

/ cos (h- 1 tE)e(E,x,x';h)xe(E)EdE(x')- k < C k , n % n {t)~ k 

Jo l^oo 



uniformly in < H < Hq. 

We call the reader's attention to the fact that these bounds decay rapidly with the 
angular momentum I ~ H . In other words, energies (in the original formulation 
of the Regge- Wheeler equation) of size <C £ 2 do not present any obstruction to the 
summation in I. This is a reflection of the expectation that any such obstruction 
should result from the local behavior of the potential around the maximum due 
to complex resonances. In the related context of the surfaces of revolutions, this 
corresponds to the principle that the growth of the constants C(£) in the decay 
estimates of [41] is determined by the local geometry of the manifold rather than 
its asymptotic behavior at the ends. In particular, if the surface contains a large 
trapping set (such as an equatorial section of a sphere) then the constants grow 
exponentially in £, rendering summation impossible. 

3. Energies close to the top, Mourre and Sigal-Soffer estimates 

For energies in the range e < E < 100 we establish a Mourre estimate which then 
allows us to invoke the semiclassical Sigal-Soffer type decay bounds of Section [B] 
Thus, let p := —ihd x , H := p 2 + V as above, and A := \(px + xp). Note that 
the Mourre estimate is shown here to hold in a neighborhood of a trapping energy 
(namely, E — 1). For notational convenience, we shift the location of the maximum 
to x max = in this section. 

Lemma 3.1. For e > and h small, there exists a fixed constant cq > so that 

(3.1) xi(H)^[H,A} X i(H)>c n X i(H) 
where xi is the indicator of I := [e/2, 100]. 

Proof. We split I = I Uh where I := [e/2, 1 - e/2] and h ■= [1 - 2e, 100]. We 
start with the latter, and write I instead of I\ for simplicity. First, 

(3.2) l - [H, A] = 2p 2 - xV'{x\ h)>p 2 - xV 
Hence, with gi being a smooth cutoff function adapted to /, 

\A 

>9i(H)( P 2 -xV') 9l (H) 

> gi{H){{p 2 - xV')F 2 + F 2 (p 2 - xV) + (p 2 - xV')F 2 + F 2 {p 2 - xV'))gi{H) 

> 9I (H) [2F(p 2 - xV')F + [F, [F,p 2 ]} + F 2 Hg I (H) + Hg^F 2 
+ 2F 2 (-xV -V)] 9I (H) 

Here 1 = F + F is a smooth partition of unity with F{x) — 1 on [— xi,xi] where 
x\ > will be a large number depending only on V. Moreover, gi is another 



gi (H)^[H,A] 9l (H) 
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function adapted to I with gigi = g%. By (|1.8|) F(— xV')F > cx 2 F 2 for some c > 
depending on x\ and the Heisenberg uncertainty principle implies that 

F(p 2 - xV')F > cF(p 2 + x 2 )F > 2c hF 2 

The uncertainty principle here is being used in the form 

Mill + Mil >2|WI|aWI| a >ftW5 
which immediately follows from the fact that \p, x] = —%H, see for example [21] . 
Furthermore, 

[F{Fp 2 ]} = -2h 2 (F') 2 

and 

gi(F 2 H~ gi + H 9l F 2 ) gi 

= 9I (F 2 (H - l)gj + (H- l)~gjF 2 + {F 2 ~ 9l + gi F 2 )) 9l 

= gi (2F(H - 1) 9I F + 2F 9I F + [F , [F, (H - 1)^]] + [F , [F, 9l ]]) 9l 

> 9I {2{1 - e)F 2 - 0{h 2 )) 9l > gi F 2 gi - Ch 2 gi F 2 gi 

where we used that || [F, [F, g]] || < ti 2 , see Lemma [3.21 below. Finally, from the 
shape of our potential V(x; H) one verifies easily that xi can be chosen such that 
—xV — V > for all |x| > X\ whence F 2 (— xV' — V) > 0. In view of the preceding, 

gi (H) l -[H.A] 9l {H)>c h g j(H) 

h 

as desired. Finally, on the interval Iq one can use ([3.21) directly since one has 
a classical nontrapping condition on energies in that range. This then gives the 
desired Mourrc estimate in that range of energies, see Theorem 1 in [50]. □ 

The following commutator bound was used in the previous proof. 

Lemma 3.2. Let F and g be smooth and compactly supported. Then 

\\[F(x),[F(x), g (H)}}\\<Ch 2 

where C = C(F, g ). 

Proof. For simplicity we show that ||[i r (x),5(i?)]|| < CH, the double commutator 
being an obvious variation thereof. By the commutator expansion formula (j6.4[) 
one has 

\\[F(x), 9 (H)]\\<C(F)\\[x, g (H)]\\ 
Now g (H) = g (H) where H := H(H + and g is again smooth and compactly 
supported. Hence one can expand with the bounded H to conclude that 

\\[F(x),g(H)}\\ < C(F,g)\\[x,H]\\ =C(F,g)\\(H + l)- 1 [x,H](H + l)- 1 \\ <C(F,g)H 
Here we used that [x,H] = -2ihp and \\p(H + l^W < \\p(l + p 2 )- 1 \\ < 1. □ 

Since we are dealing with wave rather than the Schrodinger equation, we need 
to derive a Mourre estimate for y/H rather than H. However, this is an easy 
consequence of the Kato square root formula. 

Corollary 3.3. For e > and h small, there exists a fixed constant cq > so that 

Xi(H) 1 -[Vh,A] X i(H) > c H X i(H) 
where xi is the indicator of I := [e/2, 100]. 
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Proof. One uses that 

1 f°° 

H-^xi(H) = - (H + A)- x A-i dX xi(H) 
n Jo 



whence by Lemma 13.11 

Xl (H)i[VH,A]xj(H) 

= xi(H)VH^[A,H-i}VH X i(H) 

1 poo • 

= -xi(H)VH / -[A, (H + A)- 1 ]A-5 dX VHxi(H) 
7t Jo n 

i roc ■ ^ 

= -VH / (H + X)- 1 Xl (H)-[H,A}xi(H)(H + \)- 1 X-idXVH 
n Jo h 

> -Vh / (H + X)- 1 ^ h X i(H){H + A) -1 A - ' dX \fH 

Jo 

> ~c Q h Xl {H) 

and we are done. □ 

In order to apply the time-decay result from Section [SJ we need to verify the 
basic commutator assumption (|6.ip . For the definition of ad^ we refer the reader 
to that section. 

Lemma 3.4. For any smooth function g on the line with support in (0, oo) one has 

\\ a d A (g(VH))\\<C(k,g)h k 

for all k > 1. 

Proof. For the purposes of this proof, we call any smooth function g on the line 
with support in (0, oo) admissible. First, there exists another admissible function 
g with g(yH) = g(H). Second, with H = H(H + 1) _1 for any admissible g there 
exists g admissible such that g(H) = g(H). So it suffices to consider &d A (g{H)) 
with admissible g. 

As a preliminary calculation, note that 

i[H, A] = h{2H - (2V + xV')) =: H(2H + Vi) 

i[H,A] = {H + iy l i[H 1 A]{H + l)- 1 = h(H + iy 1 (2H + V 1 )(H + l)- 1 
whence ||[.ff, < C h. At the next level, 

i[i[H,A]A] = h(2i[H,A]+i[Vi,A}) = h 2 {AH + 2V X ~ xV{) 
For H we use the general identity 

[SBS, A] = SB[S, A] + S[B,A]S + [S, A]BS 

to conclude that 
(3.3) 

i[i[H,A]A] =ih{(H+l)-\2H+V 1 )[(H+l)- 1 ,A]+(H+l)-~ L [2H+V 1 ,A](H+l)-~ L 

+ [(H+1)- 1 ,A](2H + V 1 )(H + 1)- 1 } 

Inserting 

i[(H + l)-\A] = -(H + l)- 1 i[H,A](H + 
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into (|3.3|) implies that ||z[z[ij, < Cfi 2 . Continuing in this fashion implies 

\\&d A (H)\\<C(k)K k V/c>l 

Next, we transfer this estimate to && k A (g(H)) via an almost analytic extension of 
an admissible function g. This refers to a smooth function Gn{z) in the complex 
plane of compact support such that g — Gn on the real axis and with 

(3.4) |(^G*)(z)|<CW|Im*|" 

for an arbitrary but fixed positive integer N. One then has the Helffer-Sjostrand 
formula 

(3.5) g{H) = - f (d,G N )(z)(H - z)- 1 m(dz) 

T JC 

where m is the Lebesgue measure on C, see [15j Chapter 2]. The desired estimate 
now follows from 

ad^ (<?(#)) = - f {d- z G N )(z) & A k A {{H - zy^midz) 

7T Jc 

For example, for k = 1 

ad^ff - z)- 1 ) = -(H - z)- l [H, A](H - z)^ 1 

and therefore 

\\a.d\{{H - z)- 1 ^ < C\lmz\ 2 h 
Inserting this into (|3.5[) and using (|3.4[) yields 

\\[g(H),A]\\<Ch 

The cases of higher k are analogous. The larger k is, the larger N needs to be. □ 

We are now ready to state the main decay estimate for intermediate energies. 
Corollary 3.5. One has for small H and all t > 0, as well as any a > 0, 

(3.6) \\(x)- a e i *F X i(H)(x)- a h-n < C(a)(ht)- a 
Furthermore, 

(3.7) ||(x)-«e l ^x/(^)/l|oo < Ci^h-^htr^Uxrfh 

Proof. By Corollary 13.31 and the previous lemma, we conclude from Proposition l6.ll 
that for any admissible function g (as defined in the previous proof) and any a > 

(3.8) \\(A)- a e- i *G L g(H)(A)- a f\\ a < C (hty a \\f\\ 2 

To derive (|3.6I) from this estimate, we pick another admissible g so that g(H)g(H) = 
g(H) where H = H(H + l) -1 as before. Moreover, the support of g is taken to lie 
strictly within (0, 1). Then 

(x)- a e^ X i(H)(x)- a = (x)- a g(H)(A) a (A)- a e^ X i(H)(A)- a (A) a g(H)(x)~ a 
It therefore suffices to prove that 

(3.9) \\(x)- a g(H)(A) a \\<C(a) 

The logic here is that the cutoff g(H) guarantees that H = p 2 + V is bounded, 
whence also p 2 is bounded. But then p is bounded, so A should be at most as large 
as x which justifies (|3.9[) . By complex interpolation, it suffices to prove that (|3.9j) 
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holds for positive integers a. Moreover, composing with the adjoints shows that 
this is the same as 



1 



\\(xr a g(H)(Ay a g(H)(x)- a \\<C(ay 
For example, set a = 1. Then one checks that 

1 h 2 h 2 

(A)" = 1 + -z{px + xp) 2 = 1 — — + xp 2 x = 1 — - — xVx + xHx 

Since V = 0((x)~ 2 ), it suffices to bound xHx. Let G denote the almost analytic 
extension of g as in the proof of the previous lemma. Then 

(3.10) = {x)- 1 {xg{H) + [g(H),x])H(g(H)x- {g(H), x})^)- 1 

It is clear that the terms involving no commutators are bounded. For the commu- 
tators in the second line we use the Helffer-Sjostrand formula as before, viz. 

[g(H),x} = - [ d- z G{z)[{H-z)-\x]m{dz) 

d s G(z)(H - z^lx, H}(H - z)~ l m(dz) 

(3.11) =- [ d E G{z)(H - z)- 1 ^ + l)- 1 (-2ihp){H +iy 1 (H - z^midz) 

In particular, [g(H),x] is a bounded operator. Inserting this into (|3.10l) concludes 
the argument for a = 1. For a > 1 the argument is similar. We begin by expanding 
for £ > 1 an integer 

(A) 2e = (l + (xp + px) 2 /4) e 

(3.12) = const • x mi p ni x m2 p n2 ■ ■ ■ x rn °p n ° 
where the sum extends over integer rrii , rii with 

Y n % < 2 ^ <2£ 

i i 

Moreover, using the commutator [p, a;] = — ih to move powers of p through powers 
of x, the general term in (|3.12l) may be written as x k p 2k x k where k < I. Hence, we 
need to show that 

(3.13) (x)- e g(H)x k p 2k x k g(H)(x)- £ 

with < k < £ is a bounded operator. First, the operator in (|3 . 1 3|) is nonegative, 
and moreover bounded above by 

(3.14) (x)- i g(H)x k H k x k g(H)(x)- t 

since p 2 < p 2 + V = H . Note that if we can move x k across the spectral cut-offs, 
then we are done since < k < £. To accomplish this, we start from the following 
identity, which is proved by induction: for every k > 2 

k-l 

(3.15) [x k ,H] = -2in^2x k - j - 1 px j 

i=i 

and [x,H] — —2ihp. Several comments are in order: first, domain considerations 
are irrelevant due to the cutoff g(H) which is always applied. In fact, we may 



DECAY OF LINEAR WAVES ON SCHWARZSCHILD 



2:S 



use this formally and in the end justify the procedure a posteriori by obtaining a 
bound on the i 2 -operator norm. Second, the total weight in x on the right-hand 
side of (|3.15[) is k — 1. And third, in any given term x ^~^pap we can shift the 
position of p arbitrarily using the commutator [p,x] — —2ih. To proceed, one has 

[x k ,H] = (H+l)- 1 [x k ,lT\(H+l)- 1 

so that 
(3.16) 

[g(H),x k } = - [ d- z G{z)(H - z)- 1 ^ + l)- 1 [x k ,H]{H +1)- 1 {H - z)- 1 m{dz) 

71" Jc 

Inserting (13. 15)) into the right-hand side of (|3 . 1 6[) and in view of the preceding 
comments we arrive at an expression of the form 

- { d s G{z)(H - z)- 1 ^ + iy 1 x k - 1 p{H + 1)- 1 {H - z)- 1 m{dz) 

I" Jc 

If k — 1 = we are done since p(H + is bounded. Otherwise, commuting 
x k ~ 1 through (H + to the left reduces the weight by another power. In other 
words, one obtains x k ~ 2 . Because of this reduction of the degree, the process must 
terminate after at most k commutations, and we are done with the proof of the 
first estimate (|3.6[) . 

Heuristically speaking, the second bound (|3.7I) is derived from the first by means 
of the following principle, known as Bernstein's inequality: if ip € L 2 (R) satisfies 
supp(t^) C [-R, R] (with ip being the Fourier transform), then ip <= L°°(M) with the 
bound 

|M|oo<||£||l<(2#)5|M|2 

where the second inequality is obtained by Cauchy-Schwartz followed by Plancherel's 
theorem. 

To see the relevance of this, let gi{H) with gi smooth be as above. Since p 2 + 
V < 100 on the support of gi(H), one sees - again at least heuristically - that 
also d 2 < lOOfi. -2 which restricts the Fourier support to size < CftT 1 . These 
operator inequalities can be interpreted in the sense of positive operators, or via 
quadratic forms, say. Ignoring the distinction between H and the "free "case in 
which H — Hq :— p 2 , we obtain via Bernstein that 

(3-17) \\(x)- a 9i(H)f\\oo<h-i\\(x)-<*f\\ 2 

Replacing L 2 on the right-hand side costs another h~ 5 by duality, so that one loses 
hr 1 in total over the L 2 -bound, which is what (|3.7j) claims. Note that we passed 
the weight in x through gi{H) onto / which is another technical issue, next to the 
distinction between H and Hq. 

In oder to rigorously implement these ideas it is advantageous to work with 
resolvents rather than the (distorted) Fourier transform. To be specific, we write 

(3.18) (xr a e^ Xl (H)(x)- a 

= {x)- a {l + H)- 1 {x) a {x)- a e lt -^(\ + H) 2 xi{H){x)- a {x) a (\ + H)- 1 {xy a 

Note that (1 + H) 2 xi{H) is just another cut-off. Therefore, the L 2 -decay bound 
applies to 

(x)- a e^(l + H) 2 X i(H)(x)- a 
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and it suffices to prove that 

(3.19) \\{x)- a {l + H)- 1 (x) a f\\ 00 <C{a)h-^\\fh 

which by duality then implies the corresponding L 1 — > L 2 estimate and thus im- 
plies (|3.7I) . To prove (j3. 19[) we represent the Green function, i.e., the kernel of 
(1 + H)^ 1 . in the form 

(3.20) (l + H) (x,x)-H w{ ^_ ) X [x > x ']+h W{M _ } X [x < x >] 

with W denoting the Wronskian, and where ip± are the Jost solutions to 1 + H 
which arc defined uniquely by 

-S 2 V± + Vip = v± 

tp±(x) ~ e^R as x — > ±oo 

These solutions are given in terms of Volterra integral equations in the form 

/>oo 

(3.21) iP + (x) = Tj} +>0 (x)-h- 1 e^V(y)ij(y)dy 

J X 

where ip± o : = e^f and symmetrically for By the maximum principle (or 

elementary convexity arguments - recall that V > 0) one sees that ip± > on the 
line. In view of (13.21[) therefore < ip± < ip± t o and the Green function in (I3.20P 
satisfies 

< (l + H)- 1 (x,x') < Ch^e- 1 ^ 

We used here that W(ip+,ip-) > eft -1 which follows by differentiating and/or 
evaluating (|3.21[) at x = 0. In conclusion, in order to prove (13.19[) we need to show 
that the kernel 

h- 1 {x)- a e- i ^ i {x'} a 

is bounded as an operator from L 2 — > L°° with norm < Cft - '. But this follows 
from Cauchy-Schwarz and we are done. □ 

4. Large energies 

This is comparatively easier than the other two regimes of energies. Indeed, the 
energy E is so much larger than the potential that the free case becomes dominant. 
Technically speaking, we use the classical WKB ansatz without turning points. 

4.1. The WKB ansatz for large energies. We shall use the outgoing Jost so- 
lutions f+(x,E; ft) which are defined uniquely as solutions to the equations 

-H 2 f1 (x, E- ft) + Vf+(x, E; ft) = E 2 f+(x, E; ft) 

f + (x,E; ft) ~ e ±l ~* x as x -> ±oo 

A global (at least on x > 0) representation of f+(x,E;H) is given by the WKB 
ansatz 



(4.1) f+{x, E; ft) = Eh^ T + {E ^Q~i{x, E; ft)e* /o" V^W^ + ha+ ( X} E] h )) 
where Q(x, E; ft) := E 2 — V(x; ft) and 

POO 

T + (E; ft) := / (E - \J Q(y, E; ft)) dy 
Jo 
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The prefactor is a convenient normalization, and T + guarantees the correct 
asymptotics at x = +00. This representation is valid for x > 0, which is justified 
by the bounds 

(4.2) |o + (ar, E; K)\ < (x)- 3 E~ 2 VE > 100, x > 

To obtain these estimates we start from the following equation for a(x), which is 
obtained by inserting the ansatz (|4.1[) into the defining equation for / + : 

(4.3) h(i/j 2 a)' = -ip 2 V 2 (l + ha), a(oc, E; h) = d(oo, E; K) = 
where 



and 



ip{x) := Q-i(x,E;H)e^^VQ(v^K)' 



5 (Q{x)\ 2 lQ(x) 



W= 16vQfxTJ 



16VQ(x); 4Q(x) 

ieU 2 -y(x)J + 4^ 2 -y(x) 1 w ' 

using that 1/ decays at least as fast as an inverse square. The solution of (|4.3j) is 
uniquely given in terms of the Volterra integral equation 

a(x,E;H) = [°° (1 - ei /-* VQK^) v 2 (y, £; fi) (1 + ha(y, E; h)) dy 



In addition to (|4.2I) . this integral equation implies the derivative bounds 

(4.4) \d%dia + {x,E;K)\ < (x)" 3 ^' iT 2 "* VE > 100, x > 

and all k > 0,j > 0. While these statements are routine, we now give some 
indication on how they are obtained. Write 

a(x, E; h) = ao(x, E;h) + h / k(x, y; h, E)V2(y, E; h)a(y, E; h) dy 

J X 

POO 

a (x,E;h) := / k(x,y;h,E)V 2 {y,E;K)dy 



fe(x, y; h, E):=l- e* /» d« 
To see that (|4.4[) holds for ao we expand the defining integral of do as follows: 

V 2 (y,E;h) 



a (x,E;h)= / V 2 (y,E;h)dy + ih d y 



y/Q(v,E;h) 



dy 



L y/Q(x,E;h) L V ^Q(y,E;h) 



dy 



The first two terms here satisfy the bounds (I4.4[) by inspection, whereas the integral 
involving the oscillatory kernel needs to be expanded further depending on the 
number of derivatives, i.e., the size of j + k. Note that each further expansion 
improves the decay of the integrand by one power of E and y, respectively. 
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4.2. Decay estimates in the regime of large energies. The WKB consider- 
ations of Section 14.11 imply the following decay estimate. For the definition of the 
spectral measure e(E, x, x'; h) see the low energies regime. 

Lemma 4.1. Let X>wo(E) be a smooth cutoff function supported in (100, oo) and 
equal to 1 on (200, oo). Then for all t > 0, 



sup 



sup 



h- 2 ( [ cos(tr 1 tE)e(E,x,x';h)Ex>ioo(E)dE (x')~ k f(x') dx' 
Js, J a 

<Ch- 2 (ty k J(\f'(y)\ + \f(y)\)dy 

x)- k H- 1 [ ( sm(h~ 1 tE)e(E,x,x';h)x > wo(E)dE (x')- k f(x')dx' 
Jm Jo 

<Ch~ l {t)- k [ \f{y)\dy 



Moreover, the same bounds hold as weighted L 2 — > L 2 estimates, but with (■) k 2 
instead of (-)~ k ■ 

Proof. This is essentially the same as in Section 9 of [17] . The only difference being 
the factor h. However, we leave it to the reader to check that the proofs in [17] 
easily carry over to this case as well. As for the L 2 — > L 2 bounds, they follow from 
the L 1 — > L°° ones by means of Holder's inequality. □ 



5. The proof of Theorems 11.11 and 11.21 

We begin by reducing general data to those of fixed angular momentum. Thus 

00 

(5.1) VoM = E E (M2v),^(-)>s^M 

f=o -i<j<l 

where {Xi,j}j—-e, is the usual orthonormal basis of spherical harmonics in the space 

of Y € C°°(S 2 ) with -A S 2Y = £(£ + l)Y. One has ||Y^||oo < C{£)i where 
C is an absolute constant. Now let Y be a normalized spherical harmonic with 
-A S 2Y = l(t + l)Y, and set h = t' 1 . Consider data ^[0] = (f,g)Y = (ip Q ,i^i). 
Let ip(t) denote the evolution of tp[0] under the wave equation (jl.ll) . as given by (|1.9[l 
and (II . lOjt . Then by Lemma 14. 1[ Corollary 13.51 and Proposition 12.81 provided £ is 
large, one obtains 

\\(x)' k ~^~tp(t)\\ L 2 {R . L 2 {s2)) 

< (t^h-^Ux^+^+ih^d^Jo^- 1 ^,^^^^ 

for any < k <C I and t > 0. Starting from general data ip[0] = (ipo, -0i), performing 
a decomposition as in (I5.1[) we may sum up the L 2 -bound over I 3> fc, whereas for 
the finitely many remaining £ we invoke the decay estimates from [TB] (for £ = 0) 
and [17] (for £ > 0). In this way one obtains (|1.2|) . The reason why (x)~2- 
weights are required stems from the fact the corresponding L 1 — > L°° bounds in |16] 
and [T7] need (x)~ 4 for t~ 3 decay, and then we lose another (x)~2~ due to Holder's 
inequality. On another technical note, the weights (x)~ k ~2~ for k — 3 (as in 
our case) essentially retain the orthogonality properties of the spherical harmonics 
which allows one to sum up the fixed £ bounds without any losses in £. 
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For the pointwise bounds we write ([5.1J in the form 

oo 

ip[0](x,w) = (ip ,ipi)(x,u>) = E (fe,j( x )^9ej(x))Y ej (uj) 

1=0 -t<j<t 

The evolution of these data is given by 

oo 

(=0 -t<j<t 

where iptj is the evolution of (fe,j,ge,j) under (|1.5|> . Therefore, setting fu = for 
ease of notation, and using the bound ||Y«,j||oo ^ ^ yields 

\\{x)- i i>(t)\\ L oo(ji L . L co { s2 )) 

OO 

^E^ 1 E ll^)"V/j(*)IUr 

£=0 -£<j<l 

OO 

<<r 3 E(^ E iK*>Vi(*)ik 

^=o -£<j<£ 

oo 

<(r 3 EW _i E ||<«> 4 <(-A lS )*^(i, W ),F/ )J -(a;)) sa |U i 

where we invoked the pointwise bounds of Lemma 14. 11 Corollary I3.5[ and Proposi- 
tion 12.81 for large I, and [TO] and [T7j for the remaining I. This bound can now be 
summed since 

oo 

1=0 -£<j<£ 

<\\(x)H'^s)ho(x,u;)\\ L ^ 

This implies the estimate (jl.3[) and Theorem 11.11 is proved. The proof of Theo- 
rem [O] is analogous. 



6. Semiclassical Sigal-Soffer propagation estimates 

In this section we present a semiclassical version of the abstract theory from 26] . 
Our arguments are very close to [26], but some care is required in keeping track 
of powers of h. The main result is as follows. In this section H and A are self- 
adjoint operators on a Hilbert space. H = H(H) and A — A(h) depend on a small 
parameter h £ (0, So] but with domains independent of fi. We assume the bounds 

(6.1) II ad^ (<?(#)) || <C(k,g)h k 

for all k > 0, h G (0, Hq) and smooth, compactly supported functions g on the line. 
As usual, ad A (g(H)) are the fc-fold iterated commutators defined inductively as 
<id A (g(H)) = [g{H),A} and 



a d k A (g(H)) = {ad k A 1 (g(H)),A] V k > 2. 



28 ROLAND DONNINGER, WILHELM SCHLAG, AND AVY SOFFER 

Proposition 6.1. Suppose I Cl is a compact interval so tha^ 
(6.2) xi(H)^[H,A]xi(H)>9h X i(H) 

for some 9 > 0. Both I and 9 are independent of h. Then for any smooth gi with 
support in I one has for all t G K 

(6-3) \\(A)- a e-^ gi (H)(A)- a f\\ < C (ht)- a \\f\\ 

for any a > where C depends on a, 9, gi and I, but not on h. Moreover, Hq 
needs to be taken sufficiently small depending on these parameters. 

The proof requires some preparatory work. First, recall the commutator expan- 
sion formula going back to [42], [43], and subsequently refined in [25], [45], [27) . 



(6.4) [g(H), f(A)] = J2 ^^-ad k A (g(H)) + R n 



k=l 



where /, g are smooth, compactly supported functions on the line and the error R n 
satisfies the bound 

n+2 „ 

(6.5) HiM<C»||adS(s(ff))||X; Uxf-o-^fWWdx 

k=0 J 

with a constant C n depending only on n > 1 . This error bound is obtained by means 
of the Helffer-Sjostrand formula involving almost analytic extensions of /, see [TSl 
Chapter 2]. For the expansion (|6.4[) and the error bound (|6.5p see Appendix B 
in [25], in particular (B.8) and (B.14). 

In particular, if / is of order at most p meaning that / is a smooth function on 
the line obeying the bound 



l/ (fc) (x)| < c k { x y- 



for each k > 0, then (|6.4|) can be applied with n > p. 

We now proceed as in [26] . Throughout this section, the assumptions of Propo- 
sition [nH] will be in force. 

Lemma 6.2. Let f > be or order p < 4, nonincreasing and with f(x) = for 
x > 0. Furthermore, assume that f = f£ and —fif'2 — fi where fa, f% are smooth. 
Let 1 < s < 00, a £ K, A s := (hs)^ 1 (A — a), and fix e € (0, 1] as well as n > 2. 
Then with gj as above 

(6-6) 

gi (H)i[H,f(A s )] gi (H) < s- l h9 gi {H)f\A s ) gi (H) + s" 1 ^ g I (H)f 1 (A s )g I (H) 

+ s-^ 2n - 1 -^gj(H) 

uniformly in a € M and h G (0, ho]. Here, f\ is of the same order p < 4 as f , and 
vanishes on x>0, and it depends only on f,gi and n. 



12 It would be perhaps more natural to expect Xl(H)jr[H,A]xi(H) > 9xi(H), see [20], [24] . 
The loss of an h in the lower bound is due to the fact that we establish the Mourre estimate at 
an energy which is trapping, namely the top of the potential barrier. 
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Proof. We replace H with H b := Hb(H) where b is a smooth cutoff function with 
bg = g (for simplicity, we write g instead of gi). Then 

B k :=ih- k &d k A {H b ), k>\ 

satisfy the bounds < C k by assumption for all k > 1. We begin by showing 

that 

(6.7) i[H b ,f(A s )} ~ - 8 - l (-f'(A a ))tB 1 {-f'(A.))l 

where ~ throughout this proof will mean equality up to addition of a quadratic 
form remainder rem = rem(s) satisfying the bound 

±rem(s) < a -(1+e) /i(il,) + s^ 2 "" 1 ^^ 

uniformly in h, a and with fx as above. Clearly, any term of the form rem is 
admissible for the lemma and can be ignored. Write / = F 2 and expand by means 
of (EH) 

i[H b ,f(A s )] = i[H b ,F(A s )]F(A s )+F(A s )i[H b ,F(A s )\ 

(6.8) = J2 ^- k (F {k) (A s )B k F(A s ) + F(A s )B* k F^(A s )) 

k=l 

(6.9) +s- n (R n F(A s ) + F{A s )R* n ) 

From (|6.5p and since n > 2 and F is of order < 2, one concludes that R n is bounded 
uniformly in s, a, h. We now now claim that only the term k = 1 is significant, i.e., 

(6.10) i[H b , f(A s )] ~ S -\F'(A S )B 1 F(A S ) + F(A s )BxF'(A s )) 

Indeed, we first check that the terms in (|6.8p for k > 2 are subsumed in the fx 
expression of rem. To see this, we note that 

(6.11) \(F^(A s )B k F(A s )^^)\ < \\B k \\\\F(A s )i;\\ 2 \\F^(A s )^\\ 2 

< {F(A s f^, V>) {FW (A s yi>, j>) < (fxi/j, t/;) 

provided fx is an upper envelope for both F 2 and (F^) 2 with some multiplicative 
constant. Second, for (|6.9p one uses that 

±(P*Q + Q*P) < P*P + Q*Q 

with Q := ||i?„||.s-^ 1+£ )F(A s ), P* := \\R n \\- 1 S - n +^ 1 + £ '> R n . The Q*Q expression 
is again subsumed into fx, whereas for P*P we obtain 

\\R n \\- 2 s- 2n+1+£ R n R* n < 5-2"+!+^ 

This establishes our claim (|6.10p . By assumption we can write F — u 2 , —F' = v 2 
with u,v of order < 1 whence || [Bx, u(A s )}\\ < Cs^ 1 and || [Bx, v(A s )]\\ < Cs -1 . 
Therefore, the right-hand side of (16. 10)) is of the form 

s- 1 (F'(A s )B 1 F(A s ) + F(A s )BxF'(A s )) 

= - S -\v 2 {A s )Bxu 2 {A s ) + u 2 (A s )Bxv 2 (A s )) 

~ -2s~ 1 uv(A s )B 1 uv(A s ) 

whence (|6.7|) since /' = — 2(uv) 2 . The remainder that arises here is of the /i-form 
as can be seen by arguing as in (|6.11|) . 
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To invoke the Mourre estimate (|6.2|) . we choose G smooth and compactly sup- 
ported in / and with bgG — gG — g. Then 

G(H)B 1 G(H) = G(H)^[H,A]G{H) > 9hG 2 {H) 

We now claim that 

(6.12) s~ 1 G(H)r]B 1 riG(H) ~ S - x r)G{H)B x G{H)r) 

where rj(A s ) := (—f'(A s ))^. It is clear that this claim will finish the proof. One 
has 

s^ 1 G(H)T]B 1 r]G(H) - s~ 1 r)G(H)BiG(H)rj 
= s~ 1 (riGBi [n, G] + [G,r]}BxGr]+[G,rj}Bi[r],G]) 
and there is the expansion 

n — l _fe 

[G,rj\ = £ W*^) ^ k ^A(G(H)) + s~ n R 

k=l 

The expansion of [77, G] is the adjoint of this one. To prove (|6.12[) . we observe that 
h~ k a,d A (G(H)) is uniformly bounded in k, and we also gain s _fc_1 < s~ 2 with k > 1 
since (|6.12|) is of order s _1 , and each step in the expansion gains another s _1 . The 
other issues, such as the domination by /1 etc. are very similar to what we have done 
before, and we skip them. These details are identical to those in [551 Lemma 2.1], 
see in particular the paragraph leading up to (2.10) in that reference. □ 

The following is the semiclassical analogue in this context of the key propagation 
estimate of Theorem 1.1 in 26 . 

Lemma 6.3. Let < 9' < 6 and gj be as in Proposition \ 6.1\ Let be the 

indicator functions o/R , respectively. Then for any t>0 and any m > 0, 

\\ x -(A-a- h6't)e^ gi (H) X +{A - a)\\ < C m (ht)~ m 
uniformly in h,a where C m only depends on m, 0',9, and gi. 
Proof. Define for any s > 1 

A s j := (hs)' 1 {A - a - H9t) 

Choose F > smooth, nonincreasing of order and F(x) = for x > 0. We shall 
prove the estimate 

(6.13) \\F(A^)e^ gi (H)x + (A - a)\\ < C (hty m 

To see that this implies the lemma, set s = t in (|6.13|) and note that if F = 1 on 
(—00, —6] with S > small, then 

F(A tit )x~(A-a-M't) =x~{A- a-W't), 9' := 9 - S 

Define 

4> s {t) := gi (H)f(A s>t ) gi (H), f = F 2 

ipt ■= e~' ll ^x + (A~ a)4> 

where <f> is an arbitrary unit vector. Then (|6.13|) will follow from the claim: for 
every positive integer m, 

(6.14) (Mt))t := (Vt,<M*M) < Ch- m s- m 
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uniformly in H, a and < t < s, 1 < s. Note that (<j> s (t)}t > by construction. Fix 
some to. Differentiating yields 

d t (Mt))t = (ip t ,D t <t> s (t)Tpt) = (D t <f> a )t 

D t <p s {t)= l -[H^ s {t)]+d t <j )s {t) 
n 

= g I {H) l -[HJ(A s<t )]g I {H)-s- x 6g I {H)f\A s>t )g I {H) 
By (|6.4I) . for any n > 1, 

0< (^(0)) <C n s- 2n 

The point here is that f ( - k \A sfi )x + (A - a) = for all fc > so that only the 
remainder in the commutator expansion contributes. Next, apply Lemma 16.21 with 
e = 1 to conclude that 

(6.15) D t Mt) < ^ [s-V(if)/i(^ s , t )<7/(£0 + s-^-^/Off) 2 ] 

where f\ satisfies the same hypotheses as /; in particular, it is of order zero (one can 
choose F above so that the properties of / = F 2 required by Lemma [631 are valid). 
Moreover, we fixed n much larger than to. Integrating this bound in < t < s 
therefore implies that 

(6.16) (Mt))t < C(s- 2 " + TT 1 *- 1 ) 

which implies (|6.14|) with to = 1. The idea is now to bootstrap using (|6.15|) . Indeed, 
we can apply (|6. 16|) to /i to conclude that (|6.14j) holds with m = 2. Iterating this 
procedure concludes the proof. □ 

Proof of Provosition \6.1[ This follows from Lemma IB~51 as follows. Let t > 0. First, 
write 

(A)- a = (A)- a x +(A + het) + (A)- a x -(A + het) 

The second term satisfies 

\\(A)- a x -(A + ifi0i)|| <ch- a t- a 
in the sense of operator norms. The first term we subject to the evolution: with 

3. 



a = - \Mt 



e-^ 9l (H) x + (A - a )= X -{A-a- ~-het)e^- gi{H) x + (A - a) 



X +(A - a - - A het)e-^ 9I {H) X +{A - a) 



The second term here satisfies 



\\{A)- a X + (A - a - ^h£t)e'^ 9l (H) X + (A - a)\\ < Ctr a t- a 

while the first satisfies the same bound without the weights (A)~ a by Lemma IB~3l 
which concludes the proof for positive times. For negative times one passes to the 
adjoints. □ 
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